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1. Introduction. Models of QCD must confront nonperturbative phenom- 
ena such as confinement, dynamical chiral symmetry breaking (DCSB) and 
the formation of bound states. In addition, a unified approach should describe 
the deconfinement and chiral symmetry restoring phase transition exhibited by 
strongly-interacting matter under extreme conditions of temperature and den- 
sity. Nonperturbative Dyson- Schwinger equation (DSE) models |T|, ^ provide 
insight into a wide range of zero temperature hadronic phenomena; e.g., non- 
hadronic electroweak interactions of light- and heavy- mesons 0, and diverse 
meson-meson [|J and meson-nucleon || form factors. This is the foundation 
for their application at nonzero- (T, fi) 0,0-0. Herein we describe the cal- 
culation of the deconfinement and chiral symmetry restoring phase boundary, 
and the medium dependence of p-meson properties. We also introduce an ex- 
tension to describe the time-evolution in the plasma of the quark's scalar and 
vector self energies based on a Vlasov equation. 

2. Dyson-Schwinger Equation at Nonzero-(T, fi). The dressed-quark 
DSE at nonzero- (T, /j,) is 

S~ l {Pk) = ij -pA(p k ) + ry 4 ^fc + C(p k ) + B(p k ) = ij ■ p + i^ 4 uj k+ + S(p fc ) , (1) 

where p k = (p, w k+ ), oj k+ = uj k +ifi, and u k = (2fc+l)7rT is the quark's Matsub- 
ara frequency. The complex- valued scalar functions: A(p, ui k+ ), B(p,u k+ ) and 
C(p,u k+ ), depend only on (\p\ 2 ,ul + ). With a given dressed-gluon propagator 
the solutions are determined by 

tn(~\ n~2 4 f A" nr ~x Pk-q k C{q k ) 

^)-^ = 3 J ^^-^PfeK^fe) ' (3) 

where herein we only consider models where A(p) = C(p). It is the interplay 
between the functions B and C that leads to confinement, realised via the 
absence of a Lehmann representation for the dressed-quark 2-point function |], 
0,0-0. B ^ in the chiral limit signals DCSB. 



Figure 1: Chiral order parameter: B(pk) obtained from Eq. (g), as a function 
of (T,fi): (T c = 158 MeV,/i = 0), (T = 0, fi c = 275 MeV). As described in the 
text, the phase boundary is fixed by the condition B(T,{x) = 0. 



To provide an illustrative solution of the quark DSE we employ an Ansatz 
for the scalar function characterising the dressed-gluon propagator |10| : 



D(p) = 3^5 k0 6 3 (p). (4) 

The infrared enhancement in this choice ensures quark confinement and DCSB. 
As an infrared dominant model, Eq. (Q) does not represent well the interaction 
away from {pk — Qk) 2 — and that introduces some mo del- dependent artefacts. 
However, they are easily identified and the model yields qualitatively reliable 
results, preserving features of more sophisticated studies. 

Using Eq. (H) in Eqs. (0-13) we obtain a system with two phases. The 



Nambu-Goldstone (NG) phase is characterised by dynamically broken chiral 
symmetry and confined dressed- quarks. The alternative Wigner-Weyl (WW) 
solution describes a phase of the model with restored chiral symmetry and 
deconfinement. In studying the phase transition one must consider the rela- 
tive stability of the confined and deconfined phases, which is measured by the 
(T, /^-dependent pressure difference between the two distinct phases: B(T, p) = 
-P^ng] ~~ P[Sww\- &(T, p) > indicates the stability of the confined (Nambu- 
Goldstone) phase and hence the phase boundary is specified by that curve in 
the (T, /i)-plane for which B(T, p) = . The critical line is depicted in Fig. [l]. 
The phase transition is first order for any non-zero p and second order for 
/i = 0. The model has mean field critical exponents, which is a feature of 
the rainbow-ladder truncation O]. The study of thermodynamic properties 
shows that it is essential to keep scalar and vector self-energies as well as their 
momentum dependence || |TD|| . 

Mesons are quark- ant iquark bound states and their masses are obtained by 
solving the Bethe-Salpeter equation [I2|. Here we focus on the vector channel 
and employing Eq. (f|) the eigenvalue equation for the bound state mass is JT3 



Re 



4 P± 

^2 nll^ 



a c {ul + - \Ml ± f) = 1 , 



(5) 



where v B ,c{Pl) = {B{Pl),C(Pl)}/\p 2 k C 2 {pi) + B 2 (p 2 )). The equation for the 



p-meson's transverse component is obtained with 



and in the chiral-limit yields Mi 



1 2 
o*7 . 



-Ml 

4 P- 



in Eq. 



independent of T and fi. This is 
the T = = fi result of Ref. fTl |. Even for nonzero current-quark mass, M p _ 
changes by less than 1% as T and \x are increased from zero toward their critical 
values. Its insensitivity is consistent with the absence of a constant mass-shift 
in the transverse polarization tensor for a gauge-boson. For the longitudinal 
component one obtains in the chiral limit: 



1 2 



A{fi' A - vr 2 T 2 ) . 



(6) 



The results for the medium-dependence of the p meson are summarised in 
Fig. 2. As in the case of the dressed-quark mass function, the response to 
increasing T and p, is anti-correlated: the p- mass decreases with increasing 
chemical potential and increases with temperature. This anti-correlation leads 
to an edge along which the T and p effects compensate and the mass remains 
unchanged up to the transition point. 

3. Nonequilibrium Application. The time evolution of the self energies 
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Figure 2: M p as a function of (T,fj). 



can be studied using Vlasov's equation 

dt S(p, x) + d p E(p, x)d x Sip, x) - d x E(p, x)d p f(p, x) = 



(7) 



Solving this equation is complicated for two reasons, (i) The energy is a 
functional of the scalar and vector self energies, which in general are nonzero 
and momentum- dependent. While the scalar self energy is small in the plasma 
phase due to chiral symmetry restoration, the vector self energy remains signif- 
icant ||10|| . (ii) The absence of a Lehmann representation for the dressed-quark 
propagator in the confined phase precludes the existence of a single particle 
distribution function, /, in this phase. Therefore a conventional kinetic the- 
ory is only reasonable in the deconfined phase. This situation is adequately 
represented in DSE models; e.g., Refs. |J describe a quark's (T, /^-evolution 
from a confined to a propagating mode, and Ref. |TIJ makes use of this evolved 
quasiparticle behaviour in calculating the plasma's thermodynamic properties. 



Therefore, approaching the phase boundary from the plasma domain we an- 
ticipate a discontinuous disappearance of the quark distribution function, /. 
As an illustration we employ an instantaneous interaction of the form 



D(p) = 3tt 2 7]5 3 {p) (8) 

to represent dynamics in the deconfined phase. In this case the Matsubara 
sum in Eq. (0) can be performed analytically and we obtain: 

^ - " (T^ro [1 - 2/faj:)1 ' (9) 
E °<^> = " (i +g >faVfa«) p ~ 2/(ftJ)1 ' (10) 



with the quasi particle energy: E*(p,x) = y (p*) 2 + M*(p,x) 2 , the renormal- 
ized momentum: p* — p (1 + Tr(p,x)), and mass: M*(p,x) = mo + T, B (p,x). 
As a test whether this simplification still yields necessary and qualitatively im- 
portant features, such as C ^ 1, B ^ m , in Fig. [3] we compare the momentum 
dependence obtained in the models specified by Eqs. (|4],|) in the vicinity of T c . 
Both functions are well reproduced and hence Eq. (|8|) can be used to model 
the persistence of non-perturbative effects in the deconfined domain. The so- 
lution of Eqs. (|7|.P-|T0|) provide the time-evolution of the quark self-energy and 
distribution function. 

As in the case of thermal equilibrium, the vector self energy plays an im- 
portant role. Neglecting S c and the momentum dependence of Y, B a simpler 
equation is obtained 

p 

d t f[p, x) + — rd x f(p, x) - m(x)d x m(x)d p f(p, x) = 0, (11) 

E[p,x) 

with m(x) the quark mass obtained as a solution of the gap equation in models 



without confinement. This equation has been widely studied; e.g. Refs. |L5 
However, we anticipate that the numerical solution of Eq. (0) will yield signif- 
icantly different results because of the presence and persistence of the vector 
self energy in the deconfined domain. 
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Figure 3: Momentum dependence of the quark scalar and vector self energies 
obtained in the model of Eq. (Q) compared with those in the model of Eq. (R). 
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